Temperature independent current deficit due to induced quantum nanowire vibrations 
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We consider electronic transport through a suspended volt age- biased nanowire. By coupling 
the tunneling current to a transverse magnetic field, vibrational modes of the wire are excited 
which influences the current-voltage characteristics of the system in novel ways through the induced 
electron- vibron interaction. From this analysis, we find that at high voltages the current through the 
suspended nanowire is offset from its non- vibrating ohmic value by an amount that is independent 
of both bias voltage and temperature. We also show that the corrections to the current from the 
back-action of the vibrating wire decay exponentially in the limit of high voltage, a result that holds 
even if the nanowire vibration modes have been driven out of thermal equilibrium. 

PACS numbers: 73.23.-b, 73.40.Gk, 73.63.Fg 
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Nanoelectromechanical systems (NEMS) are meso- 
scopic devices whose functionality depends on the pos- 
sibility to induce mechanical vibrations or displacements 
of one or several of their components^. Examples of such 
setups are numerous and include shuttling of single elec- 
trons and Cooper pairs 2 '^^, tuning of mechanical bend- 
ing vibrations of suspended nanowire o 6 ' 7 '^ and mechani- 
cally mediated superconducting and magnetic proximity 
effects ^ 10 ' 11 to name but a few. Crucial advantages of 
the downscaling implied by the acronym NEMS are the 
high (RF) vibration frequencies and the unprecedented 
sensitivity to external stimuli that can be achieved. This 
is due in turn to the low masses of these systems and to 
the strong coupling between mechanical and electrical de- 
grees of freedom at the nanometer length scale, see, e.g., 
Refs. dHHI and references therein. Also, nanoelec- 
tromechanical systems border the world of quantum me- 
chanics, which opens up the possibility to experimentally 
study quantum effects on the interaction between elec- 
trical and mechanical degrees of freedom in mesoscopic 
systemsi^^^. 

In this paper we will consider the nanoelectromechani- 
cal system studied by Shekhter et al^, who analyzed the 
conductance through a suspended voltage biased single- 
walled carbon nanotube (SWNT). In Ref. [Sj the au- 
thors showed that by applying a transverse static mag- 
netic field an effective entanglement between the elec- 
tronic and quantum nanomechanical degrees of freedom 
in the system can be achieved. This in turn results in 
a finite magnetoconductance at low temperatures due to 
the induced quantum mechanical suppression of the elas- 
tic electronic tunneling channel together with Pauli re- 
strictions on the inelastic channels due to the confined 
geometry of the system. At high temperatures however, 
no such restrictions on the inelastic transport channels 
are found, and the conductance exactly coincides with 
the transmission through the non-vibrating wire. 

Here we will consider the same system as in Ref. [2(| 
and in particular we will show that at large enough volt- 
ages there also exists a current deficit, as compared to 
the non-vibrating wire. Furthermore, this current deficit 
is shown to be independent of both the temperature and 
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FIG. 1: (Color online) Schematic outline of system under con- 
sideration. Electronic tunneling through the doubly clamped 
suspended SWNT under bias voltage V excite quantized vi- 
brations of the tube in the presence of a magnetic field, H. 
The effective multiconnectivity of the system leads to both 
a temperature- and bias voltage-independent current deficit 
(see text) as well as to a previously reported negative magne- 
toconductance^. Amplitude shown is greatly exaggerated. 



the bias voltage, making this system a good candidate 
device for detection of quantum vibrations in nanoscale 
systems. Also, we find that this reduction to the cur- 
rent is in general not dependent on strong coupling of 
the nanowire vibrational modes to the thermal bath as 
previously reported in Ref. 



21| where the current and 



conductance through a carbon nanotube containing an 
encapsulated fullerene was analyzed in the ballistic trans- 
port regime. Hence, the results presented are shown to 
be quite general also for oscillator distributions out of 
equilibrium. 

The system considered is shown Fig. [T] and comprises 
a doubly clamped carbon nanotube suspended over a 
trench of width L, subject to a transverse magnetic field, 
H. In Ref. [H it was shown that when the SWNT is 
biased by a voltage V the induced mechanical oscilla- 
tions of the tube lead to intermediate "swinging states" 
through which electrons can tunnel between the leads. 
By restricting the analysis to the fundamental mode 
(which gives the most important contribution) the au- 
thors showed that the system can be described by an 
effective Hamiltonian (Eq. (7) in Ref. |20|). 

Heff = Pleads + Hose + Ht , (1) 

which describes nonresonant charge transfer through the 
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suspended SWNT. In (I), the first term, 

Hl ea ds = ^ ta,kQ- a J&o,k i 



(2) 



describes the electrons in the leads; aj/ rfc [cty/r,fc] are cre ~ 
ation [annihilation] operators for electrons in state fc in 
the left /right lead with energy £;/ rj fc respectively. The 
second term in the Hamiltonian, 



H osc = hub^b, 



(3) 



describes the oscillating wire where b' [b] is a boson op- 
erator that creates [annihilates] one vibrational quantum 
and lu = (k/m) 1 / 2 is the frequency of the fundamental 
mode of oscillation with k the rigidity and m the effective 
mass of the wire (typically u> is of the order of 10 8 s _1 if 
L /an). The third term in ([1]) describes the interaction 
between the electrons and the oscillating wire, 
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J2Teff{k,k')al k a hk ,+h.c. (4) 
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Here, </> = ignxoLH/Qo is the dimensionless electron- 
vibron coupling strength, $o — h/e is the flux quantum, 
g is a geometric factor of order unity and Xq is the am- 
plitude of the zero-point fluctuation of the fundamental 
mode. T e f f(k,k') is the effective overlap integral between 
electrons in state k and k' in the left and right lead re- 
spectively. 

To calculate the charge transport through the system 
for the case when the density matrix is not in thermal 
equilibrium, we first consider the time rate of change of 
the total density matrix for the system, <r(i), which is 
given by the Liouville-von Neumann equation (see Refs. 
22H23ll24j for a similar analysis). To evaluate this we 



switch to the interaction picture with respect to the non- 
interacting Hamiltonian, H$ = Hi ea da + H OSCl for which 
the evolution of the density matrix is given by, 



ihd t a{t) = [H T (t),a], 



(5) 



where A(t) = e *Hat/h£ e —iH B t/h is any operator A in the 
interaction picture. Since we are interested in the en- 
ergy exchange between the electrons and the oscillat- 
ing wire we need only to know the evolution of the re- 
duced density matrix, which is found by tracing out the 
degrees of freedom of the leads, p(t) = Triads (o"(t)). 
By treating the electrons in the two leads as fermionic 
baths whose equilibrium distributions are virtually unaf- 
fected by the charge transfer to the wire we approximate 
a(t) ~ p{t)®<J leads where &i ea ds is the equilibrium density 
matrix for the leads. Evaluating Eq. (JSJ) to the lowest 
order in the tunneling probability enables us find find the 
equation of motion for the reduced density matrix in the 
Heisenberg picture, 



d t p{t) = -- 
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(6) 
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From this we derive the stationary equation for the re- 
duced density of the system under the assumption made 
in Ref. [2(| that the overlap integral T e ff(k,k') is inde- 
pendent of the momenta k and k', 



-[H osc , p] = |T e// | 2 [(J, + J 2 ) P + p{j\ + Jl) 



(7) 



In 0) X = V2<ft/ x o, x is the deflection operator of the 
oscillating wire and the operators Ji take on the form 
below, 



Ji = — I I de r de t I &re ixS: e iiloT e- lx± e~ iiloT 



x fr(e r )(l - fi(ei))e 



: ,r(e,-e r ) 



(8a) 



J 2 = - 



de r de, / dTe- lx£ e lHoT e lx:i e~ iHoT 



x/i(e/)(l-/r'(^))e lr(£ ^ ei) . (8b) 



Here, v is the density of states in the leads and //, r (e;,r) 
are the Fermi distribution for electrons in the left/right 
lead kept at chemical potential — ±eV/2 respec- 
tively. Multiplying Eq. ([7]) by the position and mo- 
mentum operator and tracing out the oscillator degree 
of freedom we find the following expression for the de- 
flection, (x), and momentum, (p), expectation values, 
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(9a) 
(9b) 



Eq. (|9a|) gives the force balance between the elastic 
force on the wire (left hand side) and the force in- 
duced by the charge transfer (right hand side). Since 
X oc HLg we interpret the latter as the Lorentz force, 
F = HILg. hence evaluating this equation will allow 
us to find the current, /, over the leads. Also, we note 
that J r j = Ji t 2 + 3\ 2 models electronic tunneling from 
the right/left lead respectively, and that the combination 
of these two terms effectively defines the current opera- 
tor of the system. To calculate the current we evaluate 
Eq. (|9aj) by dividing the operators J r ,i into their diag- 
onal and non-diagonal parts (subscripts d and n respec- 
tively), Tr(J r p) = Tr(J r ^p d ) + Tr(J^ n /5„), with respect 
to the eigenbasis of the oscillating wire. From this analy- 
sis we find that to the zeroth order in the operators x and 
p the operators J r j are diagonal and the expression for 
the force is proportional to the current, Jo, through the 
system (Eq. (8) in Ref. [13]). The higher order terms 
in x and p, corresponding to the non-diagonal parts of 
J Tt i, are collected in the current, I±, which, in the high 
bias limit, gives exponentially small corrections to total 
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current / (sec below). 

I = Io + h 

^ oo oo 



(10) 



i0(St+S)i 



n=0i=-n 



x / de[/ ; (e)(l-/ r (e-£M) 

/ r (e)(l-/i(e-*M)] 



7x = 1 e//l TV 



(11) 
(12) 



In Eq. (fTTjl . Go = 2e 2 W/h is the zero field conductance 
with W = \T e ff\ 2 v 2 the probability of electron tunneling 
and P(n) is the probability that the fundamental mode 
is in quantum state |n) with energy nhuj. 

Consider now the expression for the total current, 
Eqs. (JTTJ) and (Tl2^) . The current I can be understood 
by considering that the combinations of the Fermi func- 
tions puts restrictions on the allowed transmission chan- 
nels (through the Pauli principle) for electrons as they 
exchange energy with the vibrating wire. Furthermore, 
an analytic expression to (|11|) can be found in the limit 
of high bias voltage, eV ^> eVo (eVo oc maxjfcsT, Ifruj}). 
Integrating over the electronic energy this equation can 
be expressed as, 



Ihco + eV 

,/3(ehw-eV) _ l ~ e p(lhuj+eV) _ i 



7 °~E E p wm* 

n=0e=-n 

Ihw - eV 

X 



(13) 



which, in the limit of high bias voltage, reduces to an 
expression for the current that is independent of temper- 
ature and bias voltage, 

V - —d) 2, ) . (14) 



lim Iq = Go 



To evaluate the non-diagonal contribution to the cur- 
rent we expand the exponentials in Eq. © in power of 
Hqt and integrate over the electronic energies, 

Ji, 2 = - — - / dr S(t) 
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/3sinh(7rr//3) / 



(15) 



In Eq. (TJ5]), Y(H ,e™ x ,q) = [H , [H , [• . . , e^**]]] 
with q indicating the number of commutators to be eval- 
uated and (3 = (kgT)^ 1 . Evaluating ([12"]) with this ex- 
pansion we find that all contributions to the current I\ 
decay exponentially in the high bias limit [fleV 3> 1) as 
all correction terms will be of the form, 



I\ oc 



dr 



sin(eVY) 
sinh 2 (7rr//3) V C0S ( e ^ T ) 
oc {e^ v - l)- 1 
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-odd 
even 



PeV> !,<?> 2. (16) 
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FIG. 2: (Color online) Current as a function of bias voltage 
for the suspended nanotube system (arbitrary units), show- 
ing the temperature- and voltage-independent current deficit 
as compared to its non-vibrating ohmic behavior, / = GoV 
(cyan). Note that the current reduction oc cA 2 oc H 2 (red line) 
does not extrapolate to the origin. Here, coupling strength c/>2 
(black) > 0i (blue) and temperature T2 (solid) > T\ (dashed). 



Thus we find that in the limit of high bias voltage, the 
current goes as (flT)) , which differs from the ohmic behav- 
ior for the non-vibrating wire, / = GoV , by an amount 
that is independent of both the bias voltage and temper- 
ature as shown in Fig. [2l This can be understood from 
the fact that the increase in the current due to a further 
increment in the bias voltage under the conditions when 
V ^> Vo is fulfilled is not affected by the Pauli restric- 
tions on the electron-vibron energy exchange. Neverthe- 
less, the current deficit at large voltage biases is a true 
quantum-mechanical effect on transport that originates 
from the Pauli restrictions, however, these restrictions 
only affect the tunneling probability of low energy elec- 
trons close to the Fermi level. 

In contrast to the analysis of Ref. [2l[ where the distri- 
bution function P(n) was considered to be only slightly 
out of equilibrium, our results for the current deficit sur- 
vives independently of the form of P(ri), even for highly 
excited distributions. As an example, we have analyzed 
Eqs. ([7]) and (HJ) separately to order p, an analysis which 
shows that in the high bias limit the distribution function 
is indeed far from equilibrium as, e.g., the magnitude of 
the two lowest non-zero moments are2£, 



/3eV> 1 

V&TTW4>x e 
{X) fcj V 

H (eV - 2huj<j) 2 ) 



(17a) 
(17b) 



Finally, we show the current deficit as a function of the 
bias voltage for realistic experimental parameters, Fig. [3] 
Here, the influence of the multiconnectivity of the elec- 
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FIG. 3: (Color online) Plot of current deficit as a function 
of bias voltage. Here, W = 0.1, w = 10 s s~\ 0=1 (blue), 
= 0.6 (black), T — 8mK (solid) and T = 5mK (dashed). 



tron tunneling paths in the elastic channel and the Pauli- 
principle restrictions on the available inelastic channels 
are clearly visible at low bias voltages as an increasing 
current deficit. Eventually, when the new inelastic tun- 
neling channels that are added by a further increment of 
the bias voltage are not affected by the Pauli restrictions 
— which occurs at large enough bias voltages — the cur- 



rent deficit saturates to a constant value. This constant 
value of the current deficit depends on the magnetic field 
H and the mechanical properties of the nanowire through 
the parameter cf) as illustrated in Fig. [31 

To conclude we have shown that for the system orig- 
inally considered by Shekhter et aZ.— not only will the 
conductance be altered due to the multiconnectivity of 
the electronic transport through the system, but also the 
current. In particular we find that even for vibrational 
mode distributions out of thermal equilibrium the sys- 
tem is predicted to show a current offset which is inde- 
pendent of bias voltage and temperature. This is a clear 
manifestation of quantum mechanical effects on trans- 
port not previously considered which should be exper- 
imentally observable for magnetic fields of the order of 
H ^20 T (see also discussion in Ref. .20]). Also, we have 
shown that the influence of internal damping in this non- 
resonant charge transfer process decays exponentially to 
all orders in the moments of the position and momentum, 
thus making the system considered a very good candidate 
for direct observation of quantum mechanical effects in 
mesoscopic systems. 
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